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ABSTRACT 

The study of substructure in the stellar halo of the Milky Way has made a lot of progress in 
recent years, especially with the advent of surveys like the Sloan Digital Sky Survey. Here, we 
study the newly discovered tidal tails of the Galactic globular cluster NGC 5466. By means of 
numerical simulations, we reproduce the shape, direction and surface density of the tidal tails, 
as well as the structural and kinematical properties of the present-day NGC 5466. Although 
its tails are very extended in SDSS data (> 45°), NGC 5466 is only losing mass slowly at 
the present epoch and so can survive for probably a further Hubble time. The tidal shaping 
through the Milky Way potential, especially the potential of the disc, is the dominant process 
in the slow dissolution of NGC 5466 accounting for > 60% of the mass loss over the course 
of its evolution. The morphology of the tails provides a constraint on the proper motion - the 
observationally determined proper motion has to be refined (within the stated error-margin) 
to match the location of the tidal tails. 
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1 INTRODUCTION 

Within the last few years, it has become more and more obvi- 
ous that the Milky Way stellar halo is dominated by substruc- 
ture, particularly dwarf galaxies, cloud s, and tidal tails. D ata from 
the Sloan Digital Sky Survey fSDSS: lYork et al.ll200Ql) have re- 
■ vea led abundant exa mples of streams and substructure. For exam- 
\ ple. lBelokurov et all ([2006b) used a simple colour cut ^ — r < 0.4 
I to map out the distribution of stars in SDSS Data Release 5 (DR5). 
The "Field of Streams", an RGB -composite image composed of 
magnitude slices of the stellar density of these stars, showed the 
overlap of the leading and trailing arm of the well-known Sagittar- 
ius stream and the Monoceros ring very clearly. Also prominent 
was a new stream, which did not have an identified progenitor , 
and was called the "Orphan Stream" by |Belokurov et al. (2006b). 
The observational data on the Orphan S tream (Belokurov et al. 
l2007h was used by iFellhauer et"aLl (l2007h to argue that its pro- 
genitor may be the newl y-discovered disrupting dwarf galaxy 
UMa II dZucker et al.ll2006h . 

Tidal tails have proved to be an important diagnostic of 
the Galactic potential. Especially the tails of th e dissolving 
Sagittarius dwarf gal a xy (see e.g. Ilbata. Gilmore"& Irwin 1994; 
Maiewski et"aD l2003l : iHelmil l2004l : Ijohnston. Law" & Maiewskil 



2005h . which wrap around the Milky Way, are an excellent tracer 



of the strength and shape of the potential. ' Fellhauer et al ] (l2006h 
have shown with their numerical models that the bifurcation of the 
Sagittarius stream as seen in the "Field of Streams" is composed of 
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two wraps of the tidal tails and can only be reproduced if the orbital 
precession is small, i.e. if the Milky Way dark matter halo is close 
to spherical. 

Extra-tidal extensions and onsets of tidal tails have been 
claimed around a number of Galactic globular clusters in recent 
years (see e.g. 'Mevlan, Leon & CombesI l200lh . The most spec- 
tacular and convincing discovery remai ns the long and thin tail 
from th e disrupting gl obular cluster Pal 5 lOdenkirc hen et al.l200ll : 
iRockosi et al. 2002: Odenkirchen et aD E003). The tails extend at 
least 4 kpc from the cluster in the leading and trailing direction and 
contain more mass than the remaining cluster. 

Recently, two differ ent groups ( Belokurov et al.l l2006al : 
iGrillmair & JohnsonI bood) claim to have detected tidal tails of 
various extents around the disrupting globular cluster NGC 5466. 
This is an old, metal-poor ([Fe/H]= —2.22) cluster, lying at 
Galactic coordinates / = 42?15, 6 = 73? 59. In Belok urov et all 
(|2006a), the observed tails of NGC 5466 are not as long as 
those of Pal 5, stretching about 2° or 500 pc in either direction. 
IGrillmair & JohnsonI (12006 ) reported afterwards that they found ev- 
idence for a much larger extension of the tidal tails of NGC 5466. 
They claimed that the leading arm extends over ~ 30 degrees and 
the trailing arm extends at least 15 degrees, before it leaves the area 
covered by SDSS. This finding makes the tails of NGC 5466 even 
longer, but much fainter, than the tails of Pal 5. The aim of our pa- 
per is to confront these claims with theoretical expectation, as well 
as to study the survival of the tails. 

The following data for NGC 546 6 are taken 

from various sources in the literature ( HarrisI 1 19961: 

iDinescu. Girard & van Altenal Il999l : iLehmann & Schol j 
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Figure 1. Comparison of the Galactic potentials. Left: Circular velocity 
of the ML (solid, black) and DB (dashed, red) potentials. Also plotted is 
the enclosed mass (right hand y-axis) of the ML (dotted, green) and DB 
(long-dashed, blue) models. Both potentials have the same circular velocity 
(220 kms~^) at the solar radius. Right: Gravitational potential in the disc 
plane for the ML (dashed) and DB (solid) models. Also shown is the grav- 
itational potential along the 2; -axis at the solar radius for the ML (dotted) 
and DB (long-dashed) models. Even though the potentials agree well in the 
innermost part, the ML potential is much steeper in the outer parts. 



I1997I : [ Prvor et al ] [l99il). The central surface brightness is 
21.28 magarcsec"^. The total lumino sity is My = —6 .96 mag 
and the mass-to-light ratio as given by IPrvor et al J (Il99lh is ^ 1. 
Using these values, we derive a total mass of about 5 X 10"^ M©. 
The core radius has values ranging from 6.1 to 7.6 pc; the 
half-mass radius ranges from 10.4 to 13.1 pc. The most substantial 
differences in the literature occur for the tidal radius. Here, values 
are spread between 61 and 158 pc. We use this data as constraints 
on our numerical simulations. Our aim is to access possible 
initial models of this globular cluster and analyse its stability and 
evolution in different sets of Milky Way potentials. 

In the next section, we describe the setup of our simulations 
- namely, the choice of Galactic potential models, the orbit of 
NGC 5466 and finally the initial model of the cluster itself. This 
is followed by a study of the relative importance of two-body 
relaxation and disc shocking in Section [3] justifying our use of 
particle-mesh simulations in this paper. Then, in Section |4l we 
present simulations that reproduce t he shape, extent and su rface 
density of the tidal t a ils det ected by iBelokurov et al.l (l2006ah and 
lOrillmair & JohnsonI (l2006h . The properties of the remnant and 
shown to correspond to the present-day NGC 5466. Finally, we ex- 
amine how the tidal tails change as a function of the proper motion 
and hence orbit. 



2 SETUP 

2.1 Galactic Models 

Dynamical friction does not play a significant role in the evolution 
of a low-mass star cluster. So, we are able to model the Galactic 
tidal field as a smooth and analytic background potential. For the 
Galactic potential, we use one of two standard models. The first 
(hereafter ML from (M)iamoto-Nagai -1- (L)ogarithmic halo) is a 
superposition of three components. The halo is represented by a 
spherical logarithmic potential of the form 



(r) 



1 2. 
2^0 In 
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(1) 



with vq — 256 kms ^ and d — 12 kpc (and r is the spherical 
radius). The Galactic disc is modelled by a Miyamoto-Nagai po- 



tential: 

$disc(i?,2;) 



GMa 



(2) 



with Md = 10^^ Mo, h = 6.5 kpc and c = 0.26 kpc (where R 
and z are cylindrical coordinates). The bulge is represented by a 
Hernquist potential 



^bulge(r) = 



GMb 
r + a ' 



(3) 



using Mb = 3.4 x 10^° M© and a = 0.7 kpc. 

For co mparison, we als o use the Galactic potential suggested 
bylBehnen & Binney (Il998h and hereafter denoted by DB. It con- 
sists of three disc components, namely the ISM, the thin and the 
thick disc, each of the form 



/ \ ^cl / -Rm -ti^ 1^1 



R_ 

Rd 



(4) 



With i?m = 0, Eq. ^ describes a standard double exponential disc 
with scale-length Rd, scale-height Zd and central surface-density 
Ed. For the stellar discs, Rm is set to be zero, while for the ISM- 
disc, we allow for a cen tral depression by setting Rm — 4 kpc 
toehnen & Binnevlll998h . In addition to the the disc potential, we 
use the analytic potential corresponding to two spheroidal density 
distributions for the halo and the bulge in the form 



p{R,z) 
where 



(6) 



We choose the parameters of our DB model a ccording to the best- 
fit model 4 in the paper of Dehnen & BinnevI (|l998). In Fig.[T] we 
compare the two potentials. For both, the circular velocity at the 
solar radius is 220 km s~^. However, the ML model contains more 
mass within a given radius than the DB model. 



2.2 Initial Model and Orbit for NGC 5466 

As a initial model for the star cluster, we choose a lPlumme^ (Il91lh 



sphere: 



p{r) 



3Mpi 
A^Rl, 



2 \ -5/2 



(7) 



with i?pi being the scale-length of the Plummer sphere, which is 
identical to the half-light radius, and Mpi the total mass. This is a 
fairly good representation of a star cluster, especially a young one. 
Howev e r, due to tidal shaping and internal evolution at later stages, 
a Kins' (Il966l) model usually fits the photometric data better. The 
advantage of a Plummer model is that all physical quantities are 
analytically accessible. 

The initial Plummer model has a half-light radius of 10 pc, 
an initial mass of 7 x 10^ M© and is represented by 10^ parti- 
cles. The numerical set-up of the particles is performed using the 
algorithm of Aarseth, Henon & Wielen (1974). We checked that 
our initial model is able to survive for a Hubble time by com- 
pa ring our initial configuration with the dissolution times given 
in iBaumgardt & Makin ol (l2003h (see their fig. 3). As the orbit of 
NGC 5466 is most of the time located far out in the halo, it is 
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Figure 2. Comparison between the two orbits using the observationally de- 
termined proper motion. Left side (red): flipped orbit in the DB potential. 
Right side (black): orbit in the ML potential. The orbit with the revised 
proper motions in the ML potential is shown in green with an off-set of 
+60 kpc. 
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Figure 3. Bound mass versus time. Results of the NBODY 4 simulations in 
comparison with the Superbox particle-mesh simulations. The uppermost 
curve is the bound mass in the Superbox simulation. The results of the 
NBODY 4 simulations have been shifted downwards by 0.5 and 1.0 times 
10^ M0, respectively for clarity. The middle curve (green) shows the result 
of the NBODY 4 simulation neglecting stellar evolution and lowest curve 
(blue) shows the result of the NBODY 4 run with stellar evolution. Shown 
are also the linear fitting lines to the NBODY 4 results (dashed lines). 



well r epresented by the uppermost lines in 'Baumgar dt & Makinol 
(l2003h . giving us a dissolution time of a few Hubble times. 

To determine the orbit of NGC 5466, we use the positions 
and p roper motion from the literature (lHarris|[l99^ : [Binescuetal] 



ana p i 
Il999l) 



namely: 



a = 14^ 05™ 27!3 = 211?36 

3 = +28'' 32' 04'' = 28?53 

Dq = 15.9 kpc 

/j.a cos 6 — — 4.65 =b 0.82 mas yr~"^ 

^5 = 0.80 ±0.82 mas yr"^ 

Vrad — 119.7 kms""*^. 

Its Galactocentric distance is Rgc = 16.2 kpc. We transform the 
positions and velocities into a Galactocentric Cartesian coordinate 
system and integrate a test particle back in time for 10 Gyr. This 
endpoint of the backward-integration is then the starting position of 
our initial model. Even though the two potentials are quite similar 
in the innermost parts, the orbits differ in terms of perigalacticon, 
apogalacticon and number of disc crossings. In the DB potential, 
the Galaxy is less massive in the outer parts, so the cluster can 
reach an apogalacticon of 84 kpc, while in the ML case, it only 
reaches 57 kpc. The perigalactica are 5 and 6 kpc, respectively In 
Fig. [21 the shape of the orbits in the i?,z-plane is plotted (for the 
DB model, we flipped the radial coordinate onto the negative side 
to aid visibility). 



3 JUSTIFICATION OF PARTICLE-MESH SIMULATIONS 

To simulate the evolution of th e tails of NGC 5466, we use the 
particle-mesh Superbox package iFellhauer et al.l20()Qb . A particle- 
mesh code has the great advantage that we can use millions of parti- 
cles (which represent equal-mass phase-space elements rather than 
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Figure 4. Mass-loss at one perigalacticon and disc-passage of the particle- 
mesh simulation in a more detailed time resolution. The (black) parabolic 
curve shows the distance to the centre of the MW with the first (black) 
vertical line showing the time of perigalacticon. The second (red) vertical 
line shows the time of the actual disc crossing {z = 0). The (green) curve 
together with the right ordinate shows the evolution of the bound mass. It 
is visible in this figure that the mass-loss ceases after perigalacticon and 
turns into a steep mass-loss at the disc crossing. But as shown in the in-set 
(vertical lines are the same as in the main figure) the total mass-loss at this 
combination of perigalacticon and disc passage is 6 times larger than the 
steep mass-loss caused by the actual disc passage. 
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Figure 5. Mass-loss of NGC 5466 in the DB potential (left panel), in the ML potential (middle left panel), in the ML potential with the revised proper motions 
(middle right panel) and in the halo only potential (right panel). 
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Figure 6. Contour plots of the tails (the model contours have logarithmic s pacing) . The soHd green line shows the actual orbit; the red circle the size of the 
actual tidal radius. Top left: Observations using SDSS bv lBelokurov et al.l (l2006ah . Top right: Simulation using the DB potential. Lower left: ML potential. 
Lower right: ML potential combined with the revised proper motion. The tidal tails in the ML potential are more prominent than the those in the DB potential 
due to the higher mass-loss. In all models, tails and orbit (solid green line) are almost aligned. In both the top right and bottom left cases, the very inner tails 
are closer to the Galactic Centre in the leading arm and away from the Galactic Centre in the trailing arm. This is the other way round in the observations. In 
the lower right panel, the revised choice of proper motions in Eq. ( fill is used. Now the tails are a better match to the observations (compare the tails close to 
the cluster with the dashed line, which shows the 'old' orbital path). 



single stars) and trace the faint tails very accurately. However, such 
a code is often not suitable for simulations of globular clusters, be- 
cause it neglects the internal evolution due to two-body relaxation 
completely. 

The reason why Supertax is nonetheless a valid method for 
the modellin g of NGC 5466 is understood on examining the f3 pa- 
rameter (Gne din. Lee & Ostriker 1999) : 

f3 = (8) 

^shock 

Here, treiax denotes the relaxation time-scale, which amounts to ~ 
3.9 Gyr for our initial model and to '-^ 3.4 Gyr for the present state 
of the globular cluster. Additionally, tshock denotes the disc shock 
time- scale, which is the time- scale on whic h the cluster is dest royed 
by disc shocks. Using the formula from IGnedin et al J (Il999f) . we 



have 

o 2 2 

t'shock — . disc 9 5 v,"/ 

where Pdisc is the period of the disc crossings, v is the velocity with 
which the object crosses the disc, ooh denotes the ratio of velocity 
dispersion to half-mass radius rh of the object and finally Qm is the 
acceleration perpendicular to the disc. Using our simulation data, 
we derive a disc shock time-scale of about 110 Gyr. This gives a 
/3 = 0.03 ±0.01, which holds for both Galaxy models within the 
errors. The concentration 




of our initial model and the star cluster today is in the order of unity. 
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If we now place our initial model in fig. 13 of lGnedin et aT](ll999h . 
we see that it falls in the regime where shocks are more important 
than internal evolution, but also in the regime where the star cluster 
survives for a Hubble time. Still, the location of our model is close 
to the border-line (at c ^ 1 it is /3 = 0.01) where internal evolution 
becomes dominant. It is interesting to compare NGC 5466 with the 
well- studied case of Pal 5, which has c = 0.6 and [ 3 — 10. Pal 5 
will most likely be destroyed at its next disc crossing ( Dehn en et al.l 
12004) . By contrast, NGC 5466 has a good chance of surviving even 
for the next Hubble time ! 

To demonstrate that the internal evolution has no major effect, 
we perform two N-body simulations with 10^ particles in the ML 
potential using NB0DY4 (lAarsethl ll999) and compare to the Su- 
perbox results. In the first simulation, we use an equal mass for all 
particles and neglected stellar evolution. In the second simulation, 
we adopt a mass function which is present after the initial phase of 
violent mass-loss caused by the evolution of high mass stars (first 
few tens of Myr). In practice, this means another 20 % has to be 
added to the initial mass to account for the mass-loss due to su- 
pernovae, and stellar winds, as well as the stars which become un- 
bound due to this mass-loss. For the remaining stars, stellar evolu- 
tion in NBODY 4 is switched on. Figure [3] shows, by extrapolating 
the mass-loss in the direct N-body simulations linearly, that the ad- 
ditional mass-loss due to two-body relaxation and stellar evolution 
amounts at the very most to about one-third of the total mass-loss. 
The linear extrapolation of the mass-loss in this mass regime is 
justified by appeal to the work of Baumgardt & Makino ( 2003). In 
other words, disregarding the initial mass-loss when the star clus- 
ter blows away its gaseous envelope (which depends mainly on the 
star formation efficiency) and the violent stellar evolution in the 
first few tens of Myr, the dominant cause of mass-loss during the 
long-term evolution of NGC 5466 is the tidal field of the MW. 

Having established that internal disruptive processes (e.g. two- 
body relaxation) are of minor importance, we also have to prove if 
disc shocks are the major external destruction process. We there- 
fore performed a particle-mesh simulation with a halo-only poten- 
tial. The mass-loss in this case amounts to 3 % of the initial mass 
only (see Fig. [5] right panel). This is a factor of 6 less than in the 
combined potential case. But this is not yet a genuine proof of the 
importance of disc shocks. As one can see in Figs.[3]and|5]the mass- 
loss happens in short time-intervals like a step function. In Fig. |4] 
we blow up one of these short time intervals. The first (black) verti- 
cal line shows the time of perigalacticon while the second (red) line 
shows the time when the cluster crosses the disc {z — 0). While the 
mass-loss due to the tidal field ceases after perigalacticon there is 
an additional steep mass-loss starting when the cluster passes the 
centre of the disc. But shown in the small in-set in Fig.|4]this mass- 
loss is about 1 /6 of the total mass-loss at the combined perigalacti- 
con and disc passage. The general conclusion is therefore that it is 
definitely the tidal field of the disc (the perigalacticon is well out- 
side the bulge region) which causes the major contribution of the 
mass-loss, the actual disc shock when the cluster passes through the 
centre of the disc may be not that important. This finding explains 
the rather large time- scale for disc shocks (110 Gyr) of the previous 
section. 



4 TIDAL TAIL RESULTS 

4.1 The Tail Morphology and the Proper Motion 

One of the advantages of Superbox is that it has high resolution 
sub-grids, which stay focused on the simulated objects and travel 




RA 

Figure 7. Left panels: All- sky view of the tidal tails in our simulations. The 
solid green line shows the orbit of NGC 5466, its position today is marked 
with a red circle. Right panels: Maximum surface density along the tail. The 
sky is cut into square degrees in right ascension and declination and for each 
degree in right ascension the maximum surface density of all declinations 
is given. From top to bottom we show the results of the DB, ML and ML 
with revised proper motion models. 



with them through the simulation area. This is important in study- 
ing the morphology of the tenuous and diffuse tidal tails. Within 
the innermost grid, we resolve the globular cluster at a resolution 
of 1.7 pc. The grid with medium resolution is chosen to resolve the 
tidal tails close to NGC 5466 with a resolution of 16.7 pc. 

In Fig. [5] (first two panels), we plot the bound mass of our 
models in two Galactic potentials (DB and ML). In both cases, the 
mass-loss is strong in the first 2-3 Gyr and then tends to level off 
during the later stages of the evolution. The mass-loss is mainly re- 
lated to each disc crossing near perigalacticon and the mass stays 
almost constant during the rest of the orbit. This is the major rea- 
son why in the DB potential the mass-loss over 10 Gyr of evolution 
is less than in the ML potential - the star cluster has had fewer 
disc crossings. In the ML potential, the cluster loses about 18 % 
of its initial mass, whilst the cluster in the DB potential only suf- 
fers a mass-loss of 14 % (with a fluctuation of only a few particles 
out of 1,000,000). The mass-loss of the ML simulation is in very 
good a gree ment with the re s ults found in previous studies by H enonI 
(Il96lh and lLee & Ostriked imii) . However, these mass-losses are 
only lower limits, as there will be a smaller, but not negligible, con- 
tribution from internal relaxation effects. 

In the top left panel of Fig. [6] we show the data on the tails 
of NGC 5466 reproduced from Belo kurov et al.l (l2006ah . who used 
neural networks are used to reconstruct the probability density dis- 
tribution. The contours correspond to level curves of equal neural 
network output and therefore trace the star density. The tails are 
clearly visible once the extragalactic contaminants (predominantly 
galaxy clusters) have been eliminated. The tails extend ^ 4° on the 
sky, corresponding to '-^ 1 kpc in projected length. 

In the next two panels of Fig. [6l we show the projection on 
the sky of our models in the ML and DB potentials after 10 Gyr 
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Table 1. Results of the suite of simulations investigating the relation between orbit and density of the tidal tails. The columns give the following numbers: the 
absolute value of the proper motion, the proper motion in a and 5, the peri- and apogalacticon distances, the mean density of the tails, the maximum density 
in the tails, the extent of the tails, and the final mass of the cluster in units of the initial mass. 
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Figure 8. Simulations with different proper motions. The rows show from top to bottom: the projection of the tidal tails in a and 5, the distance distribution of 
the tails vs. a and finally the surface density of the tails vs. a. The columns show first the simulation with the smallest magnitude in proper motion, our best 
fitting model and finally the model with the highest magnitude of proper motion. 



of evolution. Both models show faint tidal tails, which match the 
general shape of the contours well. However, one important fea- 
ture of the data is not reproduced - the leading and trailing tails are 
well-aligned with the proper motion vector. This contrasts with the 
data, in which the inner parts of the leading tails are slightly be- 
low the proper motion vector, whilst the inner parts of the trailing 
tails are slightly above. However, the observed proper motions are 
not well determined and have large error-bars, so one possibility 



is that the proper motion of NGC 5466 should be either larger in 
right ascension or smaller in declination than the values given in 
the literature (e.g. Dinescu et al. 1999) to match the observed mis- 
alignment. We confirm this result by running another simulation 
with slightly changed proper motions, namely 

/j^aCosS — — 4.7 masyr""^ 

IJL5 = 0.42 mas yr"\ (11) 
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This orbit gives a perigalacticon of 5.9 kpc and an apogalacticon 
of 57.5 kpc. Although the change in proper motion does not make 
a significant difference to the mass-loss rate, as shown in the right 
panel of Fig. [51 it does improve the match with the location of the 
observed tidal streams much better, including the misalignment. In 
the lower right panels of Fig. [6l we see that the inner parts of the 
leading tails are now slightly below the old proper motion vector, 
whilst the inner parts of the trailing tails are now slightly above. 

4.2 The Tail Densities and Extent 

Fig. [7] shows all-sky views of the tidal tails, together with density 
profiles obtained by counting particles. The surface density of the 
tidal tails falls off along the innermost tails very steeply and stays 
at a very low density of 20-50 M© deg~^ throughout the tails. 
These low densities are ver y hard to detect, even in surveys like 
SPSS. [Grillmair & JohnsonI (l2006 l) found long, almost linear and 
very tenuous tidal extensions to NGC 5466 using a matched fil- 
ter. Although these extensions are hard to see in the SDSS data, 
they do receive some support from the simulations presented here. 
The tails of our m odel with the revised proper motion extend over 
^ 100° on the skv. lGrillmair & JohnsonI (l200 6l) claim that the aver- 
age density of the tails is about 10-20 stars deg~^, which is also in 
good agreement with out es timate. Interestingly, at the point where 
iGrillmair & JohnsonI (l2006h start to loose track of the leading arm, 
our model is close to its apogalacticon and the tails are spread out 
much wider than is the case close to the cluster. 

Although our simulation with the revised proper motion pro- 
vides a good representation of the data, it is clearly not unique. 
In particular, it is interesting to understand the variety of tidal tail 
morphology for NGC 5466, especially as forthcoming deeper pho- 
tometry will provide stronger constraints on the modelling. Ac- 
cordingly, we perform a suite of Superbox simulations to inves- 
tigate how the choice of proper motion influences the mass-loss 
and hence the properties of the tidal tails. As a constraint, we only 
used proper moti ons which a re within the la error range of the ob- 
served value (Din escu et al.l ll999) and also requir e that the orbital 
path ne ar the cluster aligns with the tails found bv lBelokurov et al.l 
(l2006ah . i.e., have the same projected orbital path as our refined set 
of proper motion. All- sky views of selected simulations are shown 
in Fig.[8]and show significant differences in the morphology and the 
properties of the tails. Table [T] gives the parameters and the results 
for the entire suite of simulations. 

The number of degrees in right ascension over which the tail 
is detectable represents a measure for the length of the tails. The 
mean density of the tails is calculated in the following way. We 
examine one degree in right ascension a and search for the highest 
surface density in the tails for each degree in declination 5. From 
these values, we compute the average surface density over the range 
of right ascension for which the tails are present. The maximum 
density given in the table is computed from the square degree of the 
tails with the highest surface density. Effects of varying distances 
are not taken into account. Table [T] shows clearly that the closer the 
orbit is to the Galactic centre the more severe is the mass-loss and 
the higher is the density in the tails. 



4.3 The Remnant NGC 5466 

Let us consider the internal properties of our remnant cluster in a 
representative simulation. We choose the one which uses the ML 
potential and the revised proper motions, shown in the lower right 
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Figure 9. Left: Surface density distribution of our model of the 
NGC 5466 remnant. Left ordinate shows pc~^, right ordinate shows 
mag arcsec"^ using the M/L-ratio from literature of 1. The inner part is still 
well fitted by the initial Plummer profile, while the outer parts are better fit- 
ted by a steeper power-law with index —4.5. The deviation is visible in the 
region from tidal radius at perigalacticon to tidal radius at apogalacticon. 
Vertical lines denote the tidal radius now (solid) and at the last perigalac- 
ticon (dashed), while dotted lines mark the grid boundaries with changes 
in resolution. Right: Velocity dispersion profile. Green open squares denote 
the 3D velocity dispersion measured in concentric shells around the cluster 
centre, crosses are the line-of-sight velocity dispersion measured in con- 
centric rings around the cluster centre. Curves show the profile of the initial 
model. Vertical Unes are as in the left panel. 

panels of Fig. [6] For this simulation. Fig. |9] shows the surface den- 
sity and velo city dis persion profiles of the final cluster. Adopting 
the data from lHarrisI (Il996) (updated values from 2003), the cluster 
has a central surface brightness of 21.28 magarcsec"^. This is in 
good agreement with our simulation, for which the central surface 
brightness is 20.6 mag arcsec"^, especially taking into account that 
our particle-mesh code neglects internal evolution, which would 
lead to higher densities in the core. Also, the half-light radius in our 
simulation is 10 pc and corresponds well with t he observed values 
of 10.4 pc (Harris 1996) and 13.1 pc JPrvor etal. 1991) . The actual 
ti dal radius in our model is 75 pc (using the Jacobi limit as given 
in'Binnev & Tremain3l 19871) and is less than the 158 pc stated in 
[Harris ( 199^ (158 pc) or 97 pc in|Lehmann & Scholz ( 1997), but 
only slightly larger than the radius of 61 pc found by Pry or et all 
([l991.) . Note that, observationally speaking, rtidai is determined by 
fitting a King (1966) profile to the surface brightness distribution, 
which does not correspond exactly to the theoretical definition. Ac- 
cording to our simulations, the tidal radius at the last perigalacticon 
was about 38 pc. 

While the surface density in the inner parts is not much af- 
fected by the mass-loss, the central velocity dispersion is reduced 
by approximately 10 %. Also visible is a rise in the line-of-sight ve- 
locity dispersion in the outer parts, which starts already within the 
actual tidal radius. This is due to the fact that all line-of-sight mea- 
surements are contaminated by unbound stars streaming in front or 
behind the star cluster. While they do not affect the central values 
because of their low number, their effect is easily measurable in the 
outer parts where the densities of the bound stars are much lower. 



5 CONCLUSIONS 

We have presented numerical simulations of the formation and evo- 
lution of the tidal tails of the globular cluster NGC 5466. We used 
direct N-body codes to argue that the evolution of the cluster is 
dominated by external effects rather than internal relaxation, and 
then grid-based codes to trace the faint tidal tails. This novel, hy- 
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brid approach is well- suited to map out the detailed morphology of 
the low-density tails of NGC 5466. 

Naively, we might expect that a low mass cluster with ob- 
served and very lengthy tails on a disc crossing orbit would not 
be able to survive f or too much longer. However, simulations by 
iDehnen et all ( l2004 l) have already shown that the disrupting glob- 
ular cluster Pal 5 has survived for at least many Gyr in a tidally- 
dominated and out-of-equilibrium state, although Pal 5 probably 
will be destroyed at the next disc crossing. Here, we have demon- 
strated that a progenitor cluster of NGC 5466, which is quite sim- 
ilar to the present cluster, could survive substantially longer, for at 
least a few Hubble times, with its extensive but tenuous tidal tails 
gradually wrapping around the whole Galaxy. 

The evolution of NGC 5466 is mainly driven by tidal 
shocks at each perigalacticon/disc crossing combination. Although 
not entirely negligible, internal effects (two-body relaxation and 
evaporation of s tars driven by post-core collapsed processes 
(lLee&Goodmanlll995l) ) play a much less important role in the 
mass-loss. It is this property which allows us to study the tidal tails 
using grid-based codes rather than the more cumbersome direct N- 
body codes. If the observationally determined mass-to-light ratio of 
1 is correct, then the initial mass of NGC 5466 is ~ 7 x 10^ M© . 
By initial mass, we do not mean the embedded mass of the star 
cluster at its formation inside a gas-cloud. If the star formation ef- 
ficiency is low, a star cluster can lose about ^ 70 per cent of its 
initial mass in stars when the gas gets blown out by high velocity 
winds or supernovae explosions. The rapid stellar evolution of high 
mass stars then adds another extreme mass-loss of ^ 20 per cent in 
the first few tens of Myr. After this initial phase of rapid evolution, 
the cluster reaches a quasi-equilibrium. This is the starting point of 
our simulations and therefore our initial mass refers to this point in 
time. 

Our numerical simulations reproduce the observational results 
of both groups who have recently studied the tidal tails NGC 5466 
with SDSS data. Mapping out the tails close to the globular clus- 
ter, [Belokurov et al.l (l2006ai) found that the leading tail emerges 
from the side pointing towards the Galactic Centre and returns 
to the orbital path from outside, while the trailing tail emerges 
from the side opposite to the Galactic Centre and returns to the 
orbital path from within. With our simulations, we showed that the 
proper motion of the globular cluster has to be smal ler in declina- 
tion an d/or larger in right ascension than reported bv lPinescu et al.l 
(I1999I) to account for the position of the tidal tails. We propose 
a new set of proper motions, fia. cos 5 — —4.7 masyr~^/i5 = 
0.42 masyr"^ for which the tail morphology is correctly repro- 
duced. This differs from the observationally determined one by 
—0.05 and —0.38 masyr"^ respectively. These changes are within 
the error margins of the observed proper motion (±0.82 mas yr~^). 

The surface density of the tidal tails falls off along the in- 
nermost tails very steeply and stays at a very low density of 20- 
50 M© deg~^ throughout the tails. These l ow densities are very 
hard to detect, even in surveys like SDSS. iGrillmair & JohnsonI 
found long, almost linear and very tenuous tidal exten- 
sions to NGC 5466 using a matched filter approach. Their work 
is supported by the simulations in this paper, which show that the 
very long (> 100°), faint tidal tails are expected. The tails in 
our simulation have roughly the same surface density as found by 
[Grillma ir & Johnson ( 2006). 

In the future, deeper photometry, radial velocities and - thanks 
to the GAIA and SIM satellites - proper motions of individual stars 
in the tidal tails may become available. Mapping out the structure 
of the tails of globular clusters and dwarf galaxies will then provide 



powerful constraints on the Galacti c potential. This work, tog ether 
with the observational pa pers of iBelokurov et al.l (l2006al) and 
IGrillmair & JohnsonI (l2006 l). has shown that NGC 5466 is a prime 
target for such studies of cold streams. Its tidal tails, though faint, 
are the longest so far claimed for any Milky Way globular cluster. 
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